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ABSTRACT 

We report the discovery of extremely broad 21 -cm Hi absorption (FWZI ~1600 km sec"') detected with the Westerbork Synthesis 
Radio Telescope in the radio source 4C37.I1 (B2 0402+379). This object has been claimed to host a super-massive binary black 
hole (Rodriguez et al. 2006). The main features in the absorption profile are two components, separated by ~1I00 km sec"'. The 
Hi absorption in 4C37.1I is unusual because it is the first case where such broad absorption is found to be centred on the systemic 
velocity of the host galaxy and not asymmetric and blueshifted as is seen in all other galaxies with broad H i absorption. Given the 
large width of the absorption, we suggest that a possible explanation for the extreme properties of the Hi absorption is that it is the 
kinematic signature of a binary black hole. If this interpretation is correct, the combined black hole mass derived from the absorption 
profile is consistent with that derived from the luminosity of the spheroid. If the broad absorption is indeed due to a binary black hole, 
this finding confirms the importance of the gaseous component in the merging process of supermassive black holes. 
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1. Introduction 

^ Galaxy mergers and accretion are common phenomena, even in 
our nearby Universe. Early-type galaxies, and in particular those 
hosting an active galactic nucleus (AGN), are believed to form 

OO through such mergers or accretion events. Because most early- 
type galaxies are known to host a black-hole (e.g., Ferrarese & 

fSj Ford 2005), it is therefore natural to expect that the two black 
holes (BH) from the progenitor galaxies will sink to the centre 

0^ of the merger remnant where eventually they will merge. During 
this process, a binary black hole system is expected to exist in the 
^ centre of the newly created galaxy. Recent numerical simulations 

• ^ (see e.g., Mayer et al. 2007 and references therein) have empha- 
sised that this process is profoundly influenced by the presence 
H I of a gaseous component, because the orbital evolution of merg- 

- - 1 ing supermassive BHs is strongly affected by friction against this 
gaseous background. Gas is indeed commonly present in the nu- 
clear regions of AGN galaxies. This gas often shows complex 
kinematics due to the effects of the AGN on the ISM (e.g., gas 
outflows; Morganti et al. 2005, Holt et al. 2008 and refs therein). 
Even more extreme kinematics could be expected in the case of 
binary black holes, where the gas in the nuclear region may be 
strongly disturbed and kinematical signatures of this might be 
observed. 

Identifying binary black holes is not easy because the dis- 
tance between the two black holes can be quite small and high 
spatial resolution is required. Indirect evidence for spatially un- 
resolved binary black holes can be found, e.g., from the mor- 
phology of the radio emission and radio jets, or from semi- 
periodic signals in lightcurves (see Komossa 2003 for a re- 
view). However, a few cases of spatially resolved systems are 
known. These include NGC 6240 (Komossa et al. 2003) and 
3C 75 (Owen et al. 1985). One candidate recently proposed is 
the super-massive binary black hole system in the radio galaxy 



4C37.11 (B2 0402+379). This radio source is classified as a 
compact symmetric object (CSO) although it has some prop- 
erties that are unusual for a CSO. Interestingly, it is one of the 
few CSOs to posses a kiloparsec-scale radio structure. On VLBA 
scales, it has two compact, variable, flat-spectrum active nuclei, 
separated by only 7.3 pc. From one of them, two relatively sym- 
metric radio lobes with a total extent of ~40 pc emanate. Because 
of this morphology, 4C37. 1 1 has been considered a candidate bi- 
nary black-hole (Rodriguez et al. 2006). 

H I absorption was earlier detected in 4C37. 1 1 against one of 
the radio lobes by Maness et al. (2004). Here we present new H i 
observations revealing that the absorption is much broader than 
found by Maness et al., and showing two components separated 
by ~1 100 km sec"'. The characteristics of the broad absorption 
in 4C37.11 are different from other, previously detected broad 
H I absorption systems, and we discuss whether the absorption 
in 4C37.11 could indicate the existence of a binary black hole 
system in this galaxy. We also present optical spectra to investi- 
gate the kinematics of the ionised gas, and to confirm the value 
of the systemic velocity which is important for the interpretation 
of the H I absorption. 

2. The broad H i absorption 

Hi observations of 4C37.11 were performed using the 
Westerbork Synthesis Radio Telescope (WSRT) on 1 Nov 2003 
(duration 5 hours) and were repeated on 4 June 2005 (dura- 
tion 1 1 hours) to confirm the detection of the broad absorption 
profile. We used two bands of 20 MHz bandwidth each, both 
centred on 1346 MHz and with 512 channels. The data were 
reduced using the MIRIAD software (Sault et al. 1995). The 
source 4C37.11 is unresolved at the resolution of the WSRT-L 
band system (about 12 arcsec). A noise level of 0.39 mJy beam"' 
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Fig. 1. H I absorption profile of 4C37. 1 1 . The kinematics of the 
optical [N II] emission-line are indicated in color; the thick verti- 
cal colored bars indicate the central velocity of the narrow (red) 
and broad (magenta) component of [Nil], while the thin hori- 
zontal colored bars indicate the FWHM of the narrow and broad 
component. The systemic velocity (v'sys) derived from the optical 
emission-lines is also shown. 



was achieved for a velocity resolution of 36 km sec (after 
Banning smoothing). The continuum flux density of 4C37. 1 1 is 
1.43 Jy. 

The observed Hi absorption profile is shown in Fig. 1. The 
extremely broad profile is detected in both observations. The H i 
absorption profile covers about 1600 km sec"' (full width zero 
intensity - FWZI) and is dominated by two peaks: one centred on 
~16000 km sec ' and a deeper, broader one centred on ~17100 
km sec"'. Fainter absorption is also detected between these two 
peaks. Despite the relatively similar bandwidth used, our H i de- 
tection is much broader than what found by Maness et al. (2004). 
While the deep H i absorption centered around ~ 17 100 km sec"' 
is detected in their VLBI observations, the remaining absorp- 
tion does not appear in the profiles presented by them, although 
no integrated profile is presented by Maness et al. (2004), so 
a direct comparison with our data cannot be done. This means 
that for the deeper absorption we know that it is located against 
the southern radio-lobe/hot-spot and extending to about 20 pc, 
while the rest of the absorption may have been resolved out by 
the high resolution o f the VLBI data. Part of the H i absorp- 
tion could be against the more extended continuum emission 
(on scales ~500 pc) that appears in the 0.3 GHz VLBI image 
(Rodriguez et al. 2006). Another possibility is that part of the ab- 
sorption is against the kpc-scale radio lobes described by Maness 
et al. (2000). However, as we will see in the next Section, our 
optical spectra indicate that the H i absorption occurs closer to 
the AGN. As discussed for the case of Cygnus A by Conway 
& Blanco (1995), velocity widths of the Hi absorption broader 
than ~150 km sec"' suggest that the absorption is unlikely due 
to gas located at large distance from the nucleus (i.e. considering 
the limited spatial scale of the continuum background, it would 
imply an overly large velocity gradient for gas associated with 
a large-scale dust-lane). This argument appears to be even more 
valid in the case of the broad absorption detected in 4C37. 1 1 . 

At the low resolution of our WSRT observations, 4C37.1 1 is 
unresolved and, therefore, the low optical depth that we derive 
for the absorption centred on -17100 km sec"' (peak r ~ 0.005) 



is a lower limit because the VLBA data of Maness et al. show 
that the covering factor is less than 1. Indeed, Maness et al. 
(2004) derive an higher value for the peak optical depth of this 
component of r ~ 0.0179 (and a column density of 1.8 x 10^° 
cm"2 for Tspin = lOOK). 



3. The broad optical lines 

In order to compare the kinematics of the neutral and ionised gas, 
we have performed new optical observations. An optical spec- 
trum of 4C37. 1 1 was taken with the William Herschel Telescope 
(WHT) on 13 January 2004 using the ISIS long-slit spectrograph 
with the 6100 A dichroic, the R300B and R316R gratings in the 
blue and red arm respectively and the GG495 blocking filter in 
the red arm to cut out second-order blue light. This resulted in 
a wavelength coverage from about 3500 to 8400 A. A slit with 
width of 1.3 arcsec has been used, aligned along the parallactic 
angle of the observations (PA 120°), i.e. almost perpendicular to 
the PA of the radio source on VLBA-scales. The observations 
were done at an airmass of 1.04-1.08 and the seeing was approx- 
imately 1.9 arcsec. The total integration time was 1800 sec for 
both the blue and red arm. Here we present results of the red 
arm (6000 - 8400 A), because [Ol], [NII-hHo'] and [S II] are the 
only identifiable emission lines in the spectrum (see also Stickel 
et al 1993). We used the Image Reduction and Analysis Facility 
(IRAF) for a standard reduction of the data. The accuracy of the 
wavelength calibration is 0.3 A. The spectral resolution derived 
from night sky-lines is 4.2 A. 

In order to analyse the kinematics of the emission-line gas, 
we fitted the profile of the [NlI-i-Ha] emission-lines. This was 
done by using the Starlink software to create ID spectra and 
fit Gaussian profiles to the emission-lines. Figure |2] shows the 
[N II-hHq'] profile across an aperture of 1 .6 arcsec in the central 
part of the host galaxy. Following atomic physics, the Gaussian 
fit to the [N II-1-Hq'] triplet was done by constraining the width of 
the [N II] doublet components to be the same and their intensity 
ratios to be 3.0 while leaving all other parameters free. A good fit 
was only obtained by using two Gaussian components (see Fig. 
2 - top). The second component is much broader and is slightly 
redshifted with respect to the first, narrower one, in agreement 
with Rodriguez et al. (2006). This result is also confirmed by 
fitting the [S II] doublet using a similar procedure (see Fig. 2 - 
bottom). 

The [Nn-i-Ha] triplet is spatially extended out to a radius 
of at least 4 arcsec (~4 kpc), but the signal outside the central 
aperture is too weak to accurately fit a 2-component model. No 
significant evidence for rotation is present in the data (although 
we note that our spectra were not chosen to be aligned with any 
kinematical axis). 

The systemic velocity derived from the average central ve- 
locity of the narrow component of the [Nll+Ha] lines in the 
central region is Vsys - 16502 + 24 km sec"' or z = 0.055046 + 
0.000078 (using the barycentric optical velocity definition). This 
redshift is in agreement with the estimate by Xu et al. (1994) and 
Rodriguez et al. (2006). 

Figure [1] shows the width and central velocity of the narrow 
and broad emission line components in relation to the H i absorp- 
tion profile. Interestingly, the broad component has comparable 
width to the absorption profile. Thus, as in many other radio 
galaxies similar to 4C37. 1 1 (i.e. young, compact radio sources), 
the ionised and neutral gas phases appear to have similar kine- 
matical characteristics. This suggests that the absorbing Hi gas 
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Spectral line 


Narrow component 


Broad component 


Av 






^center ^width 


^center ^width 






[Nil] 


16505 ± 15 544 ± 26 


16674 ± 83 1353 ± 126 


169 ± 


84 


Htt 


16500 ± 18 398 ± 36 


16845 ± 303 1634 ± 779 


345 ± 


304 


[SII] 


16469 ±18 513 ±43 


16725 ± 88 1384 ± 253 


256 ± 


90 



Notes - All units are km sec"' and values have been corrected for instrumental resolution. Av is the difference between the central velocity of the 
narrow and the broad component. Results for [S II] need to be taken with caution (because, due to the underlying stellar continuum, the intensity 
ratio of the broad component was constrained to be the same as that of the narrow component) and are merely included to show the reliability of 
our model for fitting [NII+Ha]. 
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Fig. 2. Two-component fit to the [N IH-Ha] triplet (top) and [S II] 
doublet (bottom). Details of the fit are given in Table [T] 

is located in the central (i.e. inner kpc) region of the galaxy 
where the optical spectrum was taken. 

4. Discussion: gas kinematics and binary 
black-lioles 

4.1. The unusual nature of the gas kinematics in 4C37. 1 1 

The intriguing result presented above is the detection, in a ra- 
dio source possibly hosting one of the rare examples of a binary 
black-hole, of very broad H i absorption whose main feature is 
two absorption components separated in velocity by ~ 11 00km 
sec"'. In order to be able to interpret this absorption, we first 
discuss how the properties of this absorption compare with what 
is seen in other radio sources. 

Broad Hi absorption profiles (FWZI >1000 km sec"') have 
been detected in a handful of young or re-started powerful radio 
sources (see e.g. Morganti et al. 2005 and references therein). In 
this respect, it is not completely surprising to find such broad ab- 



sorption in 4C37.1 1, a CSO which are typically considered to be 
young radio sources where the radio jet could still be embedded 
in a relatively dense medium, perhaps left over from a merger 
event. The FWZI width of the absorption profile in 4C37.11 is 
as large as that of the broadest profiles observed by Morganti et 
al. (2005). However, the feature that makes 4C37.11 unusual is 
that the absorption profile is centred on the systemic velocity. 
The 4 radio galaxies observed by Morganti et al. that have HI 
profile widths in excess of 1000 km/s, all have profiles that are 
quite asymmetric and blueshifted relative to the systemic veloc- 
ity of the host galaxy. This difference means that the absorption 
in 4C37. 1 1 may not be connected to gas outflows, or if it is, some 
key difference exists, e.g., different geometry. 

On the other hand, symmetric H i absorption profiles centred 
on the systemic velocity of the host galaxy have been found in 
many radio sources (e.g., Conway & Blanco 1995; Morganti et 
al. 2001, Vermeulen et al. 2003). However, the typical width of 
these absorption profiles is of the order of 200 km sec"', much 
narrower than what found in 4C37.11. Such relatively narrow 
symmetric absorption is often thought to be caused by circum- 
nuclear disks on the hundred-pc scale in which the rotational 
velocity reflects the galactic potential. The fact that the absorp- 
tion in 4C37.1 1 is much broader most likely means that here it 
is not caused by such a large structure. 

It is interesting to note that the ionised gas also has a broad 
component centred close to the velocity of the narrow compo- 
nent, i.e. on the systemic velocity. This is also different from 
other compact (steep spectrum) radio sources (Holt et al. 2008 
and refs therein). In these young radio sources, the broad com- 
ponents of the ionised gas are mainly found to be blueshifted 
and associated with outflows connected with the first phases of 
the nuclear activity. 

4.2. Effects of binary BH on the gas kinematics? 

The central question is, of course, what causes the broad, two- 
component absorption in 4C37.11. Maness et al. (2004) argue 
that the Hi absorption detected by them (i.e. only the red- 
shifted component) is gas expelled by the receding radio jet. 
Considering that the gas producing the H i absorption has to be 
located in front of the continuum, this can only happen if the 
outflowing gas is accelerated and dragged by the jet cocoon (so 
that it appears in front of the jet). Because now the H i absorption 
is actually found to be symmetric around the systemic velocity, 
this would then suggest that also on the side of the main jet an 
outflow of similar amplitude is present. The fact that the source 
is classified as a compact symmetric object could indeed suggest 
that the jets are oriented close to the plane of the sky, although 
this would imply very high flow velocities to explain the large 
observed (radial) component. Moreover, Rodriguez et al. argue 
that the jets are oriented < 66° from the line-of-sight, not quite 
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in the plane of the sky. A further complication is that the lo- 
cation of the blueshifted absorption is not known and may not 
occur against the main jet. Nevertheless, given all uncertainties, 
a model explaining the two absorption components as approach- 
ing and receding outflows is probably not inconceivable. 

The discussion in Sec. 14.11 shows that the properties of 
4C37.11 are quite unusual, and perhaps this means that the 
absorption is caused by a diff'erent mechanism. Maness et al. 

(2004) and Rodriguez et al. (2006) inferred from their observa- 
tions that the most likely explanation for the radio morphology 
of 4C37.11 is the presence of a binary system with two radio- 
loud, flat spectrum AGN separated by 7.3 pc. The unique char- 
acteristics of the very broad, two-component H i absorption are 
perhaps a further indication that extreme conditions are indeed 
present in this object. It is conceivable that the absorbing gas is 
part of a structure in rapid rotation as a result of the eff'ect of 
two black holes orbiting. If so, one can use the Hi velocity to 
estimate the mass of the two black holes combined. In this case, 
we would have to assume that the blueshifted part of the absorp- 
tion is located against the (northern) lobe of 4C37.1 1 at ~20 pc 
from the nucleus, i.e. symmetrical to the redshifted part detected 
in VLBl by Manness et al. (2004). The velocity separation be- 
tween the two peaks of the H i suggest that the orbital velocity is 
about ~500 km sec"'. Using these values, we find an estimated 
black holes mass of ~10^ Mq. 

This mass compares very well with the mass of the central 
compact object derived from the luminosity of the host galaxy, 
for example using the relations presented in FeiTarese & Ford 

(2005) and by Marconi & Hunt (2003). The observed /T-band 
magnitude (2MASS) implies a black hole mass of just over 10^ 
Mq. This is much higher than estimated in a similar way by 
Rodriguez et al. (2006) using a V-band magnitude. However, the 
galactic latitude of 4C37.11 is only 10.5° and the correction for 
galactic extinction is very uncertain. For example, the galactic 
V-band extinction from Schlegel et al. (1998) is almost 3 magni- 
tudes larger than that used by Rodriguez et al.. If this higher ex- 
tinction is used, a very similar value for the combined BH mass 
is derived compared to the one found from the K magnitude. 
It appears, therefore, that if the broad absorption is interpreted 
as the kinematic signature of a binary black hole, it leads to a 
consistent mass estimate. The weakness in the argument is that 
the location of the blueshifted component in not known and this 
makes a detailed interpretation of the data uncertain. It may indi- 
cate that the blueshifted part occurs on the scale of many tens of 
pc (i.e. resolved out by the observations of Maness et al. 2004). 
Even in this case the kinematics of the gas could still be aff'ected 
by the presence of the binary BH, but the kinematics would be 
more difficult to interpret and deriving the BH mass would not 
be straightforward. 

5. Conclusions 

We have presented the detection of surprisingly broad H i ab- 
sorption, mainly consisting of two components separated in ve- 
locity by ~1100 km sec"', against the compact symmetric ob- 
ject 4C37.11. Interestingly, this radio source hosts a candidate 
supermassive binary black-hole system. The characteristics of 
the H I absorption are unusual compared to H i absorption seen 
in other radio sources. The presence of a fast outflow to explain 
the broad Hi absorption cannot be ruled out, but an alternative 
interpretation could be that it may indicate the existence of a 
binary black hole in 4C37.11. If interpreted as being the kine- 
matic signature of a binary black hole, it gives an estimate of the 
mass of the combined black holes consistent with what derived 



from the luminosity of the spheroid. However, the fact that it is 
not known where exactly the blueshifted absorption is occurring 
prevents an unambiguous interpretation. It is clear that further 
high-resolution H i observations of 4C37. 1 1 are needed to deter- 
mine what is really occuring in the centre of 4C37. 1 1 . Our data 
do show that the nucleus of 4C37.1 1 is gas rich and, if 4C37.1 1 
indeed harbours a binary black hole, it underlines the importance 
of gas in the formation of such systems. 
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